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Abstract 


We consider wake-like density fluctuations produced by cosmic strings at 
the quark-hadron transition in the early universe. We show that low mo¬ 
mentum axions which are produced through the radiation from the axionic 
string at an earlier stage, may get trapped inside these wakes due to delayed 
hadronization in these overdense regions. As the interfaces, bordering the 
wakes, collapse, the axions pick-up momentum from the walls and finally leave 
the wake regions. These axions thus can produce large scale iso-curvature 
fluctuations. We have calculated the detailed profile of these axionic density 
fluctuations and discuss its astrophysical consequences. 
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I. INTRODUCTION 

If QCD phase phase transition is of first order in nature then there are many important 
consequences in the context of early universe as well as in heavy ion collisions. In the conven¬ 
tional picture of a first order phase transition, the transition happens through homogeneous 
nucleation of hadronic bubbles in the background of QGP phases. Hadronic bubbles thus 
formed will keep expanding, coalesce and eventually complete the phase transition. One of 
the important consequences in this scenario, in the context of the early universe, is the con¬ 
centration of baryon numbers inside the shrinking QGP bubbles which eventually can even 
form quark nuggets [1] due to high concentrations of baryon number inside the shrinking 
QGP phases. It was also shown that concentration of baryon number in the QGP phase can 
produce large baryon fluctuations which can survive against various dissipative processes 
[2] upto nucleosynthesis epoch and can alter the primordial abundances [3,4]. Ignatius and 
Schwarz studied [5] the presence of density fluctuations (those arising from inflation) at 
quark-hadron transition and showed that it will lead to splitting of the region in hot and 
cold regions with cold regions converting to hadronic phase first. Baryons will then get 
trapped in the (initially) hotter regions. Estimates of sizes and separations of such density 
fluctuations were made in ref. [5] using GOBE measurements of the temperature fluctuations 
in GMBR. In our earlier work [6,7], we considered the effect of cosmic string induced density 
fluctuations on quark-hadron transition and showed that such density fluctuations at QGD 
epoch can lead to formation of extended planar regions of baryon inhomogeneity. Detailed 
profile of baryon inhomogeneities and it’s effect on standard big bang nucleosynthesis were 
discussed in ref. [7]. 

Here, we study another aspect of the QGD phase transition in the context of early uni¬ 
verse if the transition is of 1st order in nature, that is, axion inhomogeneity generation during 
quark-hadron transition in the presence of cosmic string induced density fluctuations. Axion 
was first introduced to solve strong GP problem [8] and has been extensively investigated 
primarily due to its possible role as a dark matter candidate. Recently it has been discussed 
[9] that due to the larger mass of axion in the hadronic phase compared to the QGP phase, 
axions will get trapped inside the shrinking QGP bubbles at the quark-hadron transition 
(somewhat in the similar manner as baryon trapping occurs in shrinking QGP bubbles [1]). 
The axions which are initially trapped inside the QGP phase will gradually pick up momen¬ 
tum from the walls of the shrinking bubbles and leave the QGP phase [9]. These axions 
which are left behind as the QGP bubbles collapse can form iso-curvature fluctuations. The 
minimum momentum Pmin required to leave the QGP regions is given by, 

Pmin = \l - mlg) = Am (1) 

Where, mah and are the axion masses in the hadron phase and QGP phase respec¬ 
tively. 

We study this phenomenon of axion trapping when cosmic strings are present in the 
early universe. Gosmic strings are generically produced by GUT phase transitions. Further, 
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it has been recently shown that superstring theories also lead to cosmic strings with accept¬ 
able parameters [10] (see also ref. [11]). It therefore becomes important to investigate how 
such cosmic strings can affect the universe. (It is clear by now from latest WMAP data 
that contributions of cosmic strings and other topological defects to structure formation is 
insignihcant [12]). We have shown earlier that the process of hadronization will be delayed 
in the cosmic string wakes at the quark-hadron transition [6,7]. The mechanism of shrinking 
quark phases is very different in our model compared to the conventional scenario based 
on homogeneous bubble nucleation. Geometry of the collapsing interface in our case is of 
sheet like planar structure unlike spherical in the conventional case. Thus, the iso-curvature 
fluctuations produced, as the axions leave the wake region, will spread in a sheet like regions 
in contrast to the spherical clumps as studied by Hindmarsh [9]. It will be interesting to 
study the subsequent evolution of large fluctuations with such geometrical structure and 
its astrophysical consequences. The paper is organized in the following manner. In sec. 
II, we briefly review the axion production from cosmic axionic string and it’s temperature 
dependent mass following ref. [9]. Sec. Ill discusses the nature of density fluctuations as 
expected from straight cosmic string moving through relativistic fluid. Density fluctuations 
produced by the moving ’’wiggly” string is also discussed in this section. In Sec. IV, the 
dynamics of quark-hadron transition in the presence of string wakes and the mechanism of 
axion trapping in our model have been discussed. In Sec. V, we give our result. Summary 
and possible astrophysical consequences have been discussed in the dual section VI. 


II. AXION PRODUCTION FROM COSMIC STRINGS 

The presence of the so called 6 term in the QCD Lagrangian violates CP invariance. The 
most stringent bound on this 6{< 10“®) comes from the electric dipole moment measurement 
of neutron. Such extremely small value of 6 in the strong interaction theory is called ’’strong 
CP problem”. One of the most elegant solutions of this problem was first proposed by 
Peccei and Quinn [8] by introducing an axial U{1)pq symmetry. The spontaneous breaking 
of this symmetry at some scale r]a can cause 6 to settle down at some vacuum value. At the 
QCD scale instanton effects lead to explicit breaking of this U{1)pq symmetry forcing 6 to 
settle down at the true vacuum 9 = 0. However, due to spontaneous breaking of U{1)pq 
symmetry, the phase of the held can wind around the vacuum manifold non-trivially in the 
physical space, forming the so called axionic string. After formation the strings are stuck 
in the plasma and are stretched by the Hubble expansion. During this time the dominant 
mechanism of dissipating energy is via heat due to the large frictional force exerted by the 
background plasma. However, with time the plasma becomes dilute and below a temperature 
T* which is given by [13], 

r. = 2 X (2) 

the strings move freely. The corresponding time f*, can be obtained from the following 
time-temperature relation, 

rj,2 _ 0-3 ^pl /gN 
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and is given by, 




0-3 ^Pi ^ ^ io-2i5eci 


vr*T^ 


.W^GeV, 
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(4) 


Here = 106.75 is the number of degrees of freedom relevant at temperature T* (taking 
all the species in the standard model). After the string will lose energy through radiation 
of (pseudo) Goldstone bosons which are called axions. The axion is not truly a Goldstone 
Bosons. It will acquire mass due to the axial anomaly [14] once the instanton effect turns 
on. In the high temperature phase (T > Tc), the mass of axions arises from the ^-dependent 
part of the free energy density which have been calculated by Preskil et al. [15] in the 
dilute-instanton- gas approximations. For three quark flavors appropriate for T not too 
high compared to Tc, they have determined the mass as, 


m, 


aq 


(T) ~ 2 X 10“^; 
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Where, m„, mg are the current quark masses of u, d and s quark respectively and 
A = Aqcd ~ 200MeG is the QGD scale. Taking critical temperature Tc ~ 150MeG and 
putting the values of u, d and s quark masses (see ref. [9] for details), we get the axion mass 
at the onset of QGD transition and can be written in terms of axionic string formation scale 
as. 


rriaqiTc) = 5.3 x 10 ^eH( 




Va 


( 6 ) 


At low temperature in the hadron phase, the axion mass is calculated from the current- 
algebra technique which gives the value [15], 
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Here, ttItt and are the mass and decay constant of pions respectively. The axion 
mass being temperature dependent at high temperature phase as in Eq.(5), it’s very light 
and ultra-relativistic at very high temperature. It becomes dynamically significant at the 
temperature at which compton wave length of the axions fall within the horizon. This is the 
time (t), at which axion string will not be able to oscillate sufficiently to radiate axions due to 
damping of domain walls attached to each string. From the time temperature relation Eq. (3) 
one can get the time f, from the following relation (neglecting the logarithmic dependence 
on T in Eq.(5)), 


magit) = maqiQi^f = 4(t) = (2t) 

Zr 


-1 


( 8 ) 


Where, maq{tc) is the mass of the axion at the onset of QGD phase transition obtained 
from Eq.(5). Thus, one gets the time t as. 


i=\[muq{QG] 1.6x10 '^ec(Y^^^^)i/='(^)2/3 


(9) 
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Here, we will study the fluctuations in the density of axions which are produced through 
the radiation from the axionic string during the era, f t as in ref. [9]. During this 
time one can assume that string network will enter into the scaling regime which has been 
discussed extensively in literature (see, ref. [16]). Under the scaling assumption the spectrum 
of axions in the comoving momentum range k* and k can be considered to be flat [9]. We 
will take the number density distribution of axions resulting from axionic string decay in 
the comoving momentum range k*{= g^g given in ref. [9] as, 

rikdk = R{i)Gr]lln{r]ai)pcdk/k‘^ (10) 

Where, R(t) is the scale factor at time i, and pc is the corresponding critical density of 
the universe. We will consider the above density distribution with suitable dilution factor 
due to expansion of the universe inside the wake-like overdensity density regions produced 
at QCD epoch as discussed in the following section. 


III. DENSITY FLUCTUATIONS ARISING DUE TO LONG COSMIC STRINGS 

In this section we will review the formation of density fluctuations due to cosmic string 
moving through relativistic fluid. First, we will discuss the case for a long straight string 
assuming it does not have any small scale structure. Though, simulations [17] on the string 
network in the expanding universe show that the long cosmic string possesses substantial 
amount of small scale ’’wiggly” structures running along the strings. The presence of such 
wiggles makes the string heavier and slows down the motion of the string as a whole. How¬ 
ever, as we discuss below, even the presence of such wiggles along the string will not alter 
the order of magnitude estimate of density fluctuations produced at QCD epoch. Therefore, 
for simplicity of presentation we will consider the case of straight (i.e., with no small scale 
structure) strings first and then we briefly discuss for the case of wiggly strings. 

To discuss the generation of density fluctuations by moving straight string, we take the 
metric describing the space-time around the string (lying along z-axis) as [18], 

ds'^ = dt^ — dz^ — dr‘^ — (1 — 4:GpYr‘^d'ip'^. ( 11 ) 

Here, p is the mass per unit length (equal to string tension) of the straight string. This 
metric can be put in Minskowski form by defining angle tjj' = (1 —4G/i)'^, where new angle tjj' 
varies between 0 and (1 — 4G/r)27r instead of 0 and 27r. A test particle at rest with respect to 
the string does not feel any gravitational force. However, when the string moves with velocity 
Vst-, the nearby matter gets an impulsive velocity along the direction of the surface swept by 
the moving string. For collisionless particles, the resulting impulsive velocity [16] is Vimpui = 
ATiGpVst'yst- This effect is responsible for the formation of wake [19] and for collisionless 
particles, this leads to planar structure of density fluctuation. In this case, the density 
flnctnations will be of order nnity with opening angle of the wake being eqnal to the deficit 
angle SttG/t. However, at QCD epoch the universe was dominated by the relativistic plasma. 
So, to study the density structure of the wake at that epoch, the description of matter 
should be taken in terms of a relativistic flnid rather than collisionless particles. Generation 
of density flnctnations due to a cosmic string moving throngh a relativistic flnid has been 
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analyzed in the literature [20-22], where it has been shown that for supersonic velocities 
of string, a wake-like shock is formed behind the string. We will follow ref. [22], where 
the equations of motion of a relativistic fluid are solved in the string space-time (Eq.(ll)), 
and both subsonic and supersonic flows are analyzed. One hnds that for supersonic flow, a 
shock develops behind the string, just as in the study of ref. [20,21]. Following ref. [22], the 
resulting density fluctuations can be expressed in terms of four-velocity and given by. 


5p IQt^G p,u]{l + ul) 

_ i_ 

P ?>Us^U^f — ul 


sinOw — —, 


( 12 ) 


where «/(= Vf/^1 — vj) and Us{= Vg/— v^) are the four velocity of fluid (in the string 
rest frame) and four velocity of sound respectively, with Vg being the sound speed. 6^, is 
the angle of the shock. In this case, when string velocity Vf is ultra-relativistic, then one 
can get strong over densities (of order 1) and the angle of the wake approaches the deficit 
angle ~ SnGp,. To discuss the effect of density fluctuations on QCD phase transition we 
will follow our previous work [6,7] and take the sample value of density fluctuations and the 
angle of the shock as. 


9^ ~ 20°, 5p/p ~ 3 X 10“^ (13) 

These values are obtained from Eq.(12) for string velocity of 0.9, sound velocity Vg = 
1/v^ and Gp, ~ 10“° (for 10^° GeV GUT strings). 

Here, we should mention that the axion being very weakly interactive [9], they can be 
considered to be a fluid consisting of collisionless particles. Therefore, description of wake 
formation for collisionless particles as mentioned above can be operative for axions also and 
can lead to formation of axionic density fluctuations axion of order unity with angle of 
the wake, 6^ ~ SttG/t. These over dense axions could be concentrated within very thin 
sheet like region of thickness SttG/t du ^ 1cm {dn ^ 10 km being the horizon size at 
Tc = 150 MeV). (One has to properly account for the wavelength of axions.) However, 
we will see below that the magnitude of the density fluctuations produced through axion 
trapping (which is our main topic of discussion in this paper) by the collapsing interfaces 
will be very large compared to the above density fluctuation. Essentially, this thin sheet 
like region will be trapped initially inside the wake of larger thickness which is produced 
by the formation of shock by moving string through relativistic fluid (as discussed above). 
Ultimately, axions will also escape from this thin wake and contribute in producing overall 
iso-curvature fluctuations. However, their effect will be very insignihcant. Even the presence 
of wiggles along the string as we discuss below, will not change the order of magnitude 
estimate of such axionic density fluctuations. Therefore, to make the physical picture of 
axion trapping simpler, we will not consider these kind of fluctuations in this work and we 
will focus only on iso-curvature fluctuations produced through trapping of axions inside the 
collapsing interfaces of the wake. 

So far we have discussed density fluctuations produced by the straight string assuming 
it does not have any small scale structure. But, as we mentioned above, the string has 
small-scale wiggly structures [17] running along the string. The presence of such wiggles 
increase the average energy density of the string and makes the strings to move slower. For 
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physics at large scales (compared to the scale of wiggles 1^, which is expected to be set by 
the gravitational-radiation from the string [23] and given as, ~ TGfit where F ~ 100), 
the effect of wiggles is taken into account by dehning the effective mass per unit length 
and the string tension T^. Unlike the case of straight string, these quantities for perturbed 
string are not equal and they are related by the equation of state [23], fiwTw = Due 
to the presence of these small scale structure, the static wiggly strings develop Newtonian 
gravitational held and a test particle at rest with respect to such perturbed string experiences 
attractive gravitational force. The metric around the wiggly string and density structure of 
the wake for collisionless particles have been studied and discussed in the literature [16,23]. 
It turns out that, for moving wiggly string, total impulsive velocity imparted on the particles 
towards the surface swept out by the string can be written as [23], 


^impul st "F 


27rG{n^ - Ty,) 
'^st'yst 


Where, hrst term is the usual velocity impulse arises due to the conical structure of the 
space-time around the wiggly string, with fi being replaced by the effective mass per unit 
length fiw The second term arises due to the Newtonian gravitational held developed by 
the wiggliness of the string. Using the value of ~ 1.9/i and ~ 0.5/i (see ref. [16]), 
we see from Eq.(14) that if the string moves with ultrarelativistic speed, then the hrst term 
will dominate over the gravitational term. The rms velocity of the wiggly string on the 
scale of smallest wiggles found out to be, [24] Vgt ~ 0.6. For this case, gravitational term 
can be neglected. In this case, density huctuations for collisionless particles will be order 
of unity with opening angle, 6^ — SttG/Tu,. However, the velocity of the string obtained 
by taking average over a correlation length comes out to be very small {vgt ~ 0.15) for 
which gravitational ehect should be taken into account. Formation of density structure of 
wake for collisionless particles have been studied in ref. [23] assuming the particles remain 
at a distance far from the string and the particles will experience only averaged ehect of the 
wiggles (in this case, relevant velocity is, Vst ~ 0.15 ). In the string rest frame (string lying 
along z-axis), if the particles moves with velocity Vgt in the (-l-ve) x direction, then following 
ref. [23], resulting density huctuations can be written as. 


( ^axion — 1 “1“ 4G(^.y; 

P 



x-Xo 

Xo 


(15) 


Where, Xq is of order of the inter-string separation. Note, setting the second term equal 
to zero, one recovers the result for density structure of the wake for straight string. If 
we put the value of wiggly string parameter in Eq.(15), we see that order of magnitude 
of density huctuations will not be ahected compared to the case of straight string. We 
will see later that due to the trapping mechanism, axion density can increase in the wake 
by many orders of magnitude. We thus conclude that generation of density huctuation 
by axions (being collisionless) by above phenomenon is insignihcant compared to the case 
where axionic huctuations happens through trapping mechanism inside the wake produced 
by relativistic huid. Therefore, we will ignore this ehect in our subsequent discussion. 

Now, we discuss the ehect of wiggliness on the value of density huctuations (Eq.(13)) 
produced by moving string through the relativistic huid. Study on formation of shock 
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by moving wiggly string has been done in the work of ref. [23] under some simplifying 
assumptions. In their study, authors have treated the fluid as a non-viscous, non-relativistic 
fluid. As has been argued in ref. [23], the assumptions of non-relativistic fluid can be good 
approximation for the string velocity Vgt ~ 0.2 which is relevant for the distance scale larger 
than the scale of the wiggles. Formation of wakes for wiggly strings has been discussed 
in ref. [23] after recombination. It is also mentioned in ref. [23] that similar wake should 
arise due to wiggly string even at earlier stages when string moves through the relativistic 
plasma, even with Vgt = 0.15 for the wiggly string. If we follow the approach of ref. [20-22], 
then Vst = 0.15 is subsonic and no shock can form. However, this does not represent the 
actual physical situation. This is because at smaller distance scales (compared to the scale 
of the wiggles), the straight segments of the string move with Vst — 0.6 (ref. [24]). This is 
supersonic motion and shock will form. Due to random velocities of different segments of the 
string, the individual shocks will combine to give some resultant shock and hence, wake for 
the total length of the (wiggly) string. Thus we will assume that even with the wiggles, shock 
and wake formation results. For these wakes we will use Eq.(12) with Vst = 0.6 and mass 
per unit length jj, for the straight string (as individual segments are straight). However, now 
the length scale of the shock (wake) will be governed by the size of each straight segment. 
The length of the wake (in the direction away from the string) will still be governed by the 
average inter-string separation. That will not be affected by the presence of the wiggles. 
(Recently, the evolution of wiggly string network has been studied [25] on flat space-time. 
The results suggest that even in the presence of wiggles the string network obey scaling 
solution. Authors in ref. [25] also mentioned that the average inter-string separation also 
scales with time. Though above simulation has been done on flat space-time, the author 
believes that such scaling solution is expected to exhibit in expanding universe also. If this is 
the case, then on an average inter-string separation for long wiggly string will not be altered 
from the straight string case.) However, the width of the wake (along the string direction) 
will now be truncated down to the average size of individual straight segment. Note however 
that, this simply changes the geometry of a planar wake (for a completely straight string) 
to a collection of strips with each strip corresponding to individual straight segment of the 
wiggly string. Even the geometry of these strips may not be planar due to rapidly changing 
velocity of string segments. For a correct treatment, one should properly account for the 
different velocity and orientations of these strips, however, for present we will ignore these 
complications. Following Eq.(12), the opening angle of the wake and density fluctuation for 
each straight segment of the wiggly string can now be written as (using Vgt = 0.6 and sound 
velocity Vg = 1/v^), 


6*^ ~ 70°, 5p/p~4xl0 (16) 

As the individual shock will combine to give resultant shock, we expect that the overall 
volume of the wakes will remain same (upto some factor) as for straight string case. Here, 
we should again emphasizing that the above value of density fluctuation and opening angle 
of the wake (i.e., Eq.(16)) is derived from Eq.(12), where it is assumed that the fluid flow 
to be uniform. In reality, due to the wiggliness of the string the fluid will experience the 
rapidly changing directions of the wiggles which causes acceleration of the fluid in different 
direction. So, in proper treatment of shock formation by wiggly string, one should take the 



non-uniformity nature of the fluid in the analysis. However, due to lack of such complicated 
proper analysis, we will use (mainly) Eq. (12) and Eq.(16) and discuss the effect of such 
fluctuations on QCD phase transition. 


IV. EFFECT OF STRING WAKES ON QUARK-HADRON TRANSITION AND 
AXION TRAPPING MECHANISM IN OUR MODEL 

Here, we will first briefly discuss the dynamics of QCD phase transition in the context 
of early universe. In the conventional scenario of first order phase transition, the transition 
happens through the nucleation of hadronic bubble in the QGP background. These bubbles 
then grow, coalesce and hnally convert the whole QGP phase into hadronic phase. However, 
since the critical size of the bubble beyond which it can grow is quite large at the temperature 
closed to Tc and also the nucleation rate too small at such temperatures, universe has to 
supercool to a temperature T^c to begin the nucleation process. The amount of supercooling 
ATsc depends on the surface tension a and the latent heat L. We take sample values of 
these parameters as obtained from the lattice results [26], a ~ O.OIST^^ and L ~ 3T^. As 
we discussed in ref. [7], for smaller value of surface tension our results remain applicable as 
the supercooling would be smaller (as has been argued in the literature [27], see also ref. 
[28].) and cosmic string induced density fluctuations will have more prominent effect on 
the dynamics of quark-hadron transition. Taking the values of these parameters one can 
estimate the amount of supercooling to be [5,29] 

AT,, = 1 - ^ ~ 10-^ (17) 

for critical temperature T, = 150 MeV. 

As soon as the universe supercool down to temperature T,,, hadronic bubbles will start 
nucleating in the background of quark gluon plasma phases. The duration of this nucleation 
process is very small and lasts only for a temperature range of AT„ = (1 — ^) ~ 10“®, for 
a time duration of order Atn — {tn is the Hubble time). The bubbles which have 

already been nucleated will keep expanding and it will release latent heat which will prevent 
further nucleation process. After the nucleation processes stop, the universe enters into the 
so called slow combustion phases [1]. However, in our model [6,7] this slow combustion 
phase is very different from the standard scenario. Here the phase transitions happens in 
the presence of density fluctuations (which will transform into temperature fluctuations) 
which is produced as the string passes through the relativistic plasma. In this scenario, as 
we discuss below, the hadronization inside the wake will be delayed while outside the wake 
the universe will enter into slow combustion phases. To understand this slow combustion 
phase in our model we take density fluctuation from Eq.(13), which will transform into 
temperature fluctuation of magnitude AT^ake = ST/T ~ 10“^. Since, this temperature 
fluctuation is larger than AT„, there will be no nucleation inside the wake while nucleation 
process will get completed outside the wake. Thus, the region outside the wake will enter 
into slow combustion phase, while overdensity region will till remain in QGP phase. For 
this, the overdensity in the wake should not decrease in the time scale Atn- To calculate 
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the time scale tshk of the evolution of the overdensity we take the typical average thickness 
of the wake as, 


dshk — sinOyj du — 3km. (18) 

Where, dn — 30 km is the horizon size at Tc = 150 MeV. (Note, as we have discussed 
earlier, the shock parameters of Eq.(16) apply only for a straight segment of the string. 
However, the resulting wake may still extend upto typical inter-string separation. Thus the 
overall volume covered by the wakes of different segments of wiggly string could be of same 
order as that of a wake from a straight string. We take this to be the case.) Now, typical 
time scale for the evolution of the overdensity will be governed by the sound speed Vg- We 
take Vg — 0.1 suitable at temperature closed to Tc (see ref. [5,30]). Thus, we get the time 
scale as, 

tshk ^ ^ 20 sinO^ tn (19) 

Vg 

which is large compared to the duration of the nucleation process ~ 10~^tH- It 
is also much larger than Attmsni— 34fisec) during which the quark phase is completely 
converted to the hadronic phase in the region ontside the wake. Thus, one can conclnde that 
the region ontside the wake enters the slow combustion phase before any signihcant bubble 
nucleation can take place in the wake region. Heat released by the bubbles outside the wakes 
will prevent bubble nucleation everywhere. It is then possible that the region outside the 
wakes will be converted to the hadronic phase while inside the wakes QGP phase remains. 
Further completion of the phase transition will happen when the interfaces, seperating the 
QGP phase inside the wakes from the outside hadronic regions, move inward from the wake 
bonndaries. These moving, macroscopic, interfaces may trap most of the baryon nnmber 
in the entire region of the wake (and some neighborhood) towards the inner region of the 
wake. Finally the interfaces will merge, completing the phase transition, and leading to a 
sheet of very large baryon nnmber density, extending across the horizon. Detailed prohle of 
baryon inhomogeneities prodnced throngh these collapsing interface had been estimated in 
our previous work [7]. 

Here, we want to discuss the mechanism of axions which will get trapped initially inside 
the wake and subsequently leave after acquiring required momentum form the collapsing 
walls. For that we will concentrate on a single wake like overdensity region. Typical average 
volume of such region will be governed by the total nnmber of long strings in a given horizon. 
From nnmerical simulation [24], the nnmber of long strings is expected to be about 15. (As is 
mentioned earlier, simulation [25] on wiggly string network snggests that even in the presence 
of wiggles on long string, the network obey scaling solntion. Though above simulation has 
been done on flat space-time, the author believes that such scaling solution is expected to 
exhibit in expanding nniverse also. If this is the case, then the number of long strings as we 
have taken for straight string case will remain applicable even for wiggly string also.) If the 
string wakes are reasonably parallel, then they will span most of the horizon volume, as the 
average thickness of a wake will be order of 1-2 km. In such a situation, the hadronic phase 
will first appear in the regions between the wakes, which may cover a very small fraction of 
the horizon volume initially. If / be the volume fraction occupied by the QGP phase then 
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the initial value of this fraction / will be close to 1. / will then slowly decrease to zero 
as the planar interfaces (formed by the coalescence of bubbles in the region in between the 
overdense wakes) move inward, converting the QGP region inside the wake into the hadronic 
phase. Certainly, the actual situation will be more complicated than this, with string wakes 
extending in random directions, and often even overlapping. In such a situation, even the 
initial value of /, when bubble coalescence (in the regions between the wakes) forming planar 
interfaces, may be smaller than 1. Note that, for wiggly string case, the wakes formed by 
straight segments of the string will have different orientations depending on the direction 
of each straight segment. In this situation, distribution of string wakes will be even more 
complicated, hence in determining the exact initial value of /. However, for simplicity, we 
will take the initial value of / to be almost 1 as is done in ref. [7] and focus on the region 
relevant for only one string, covering about 1/15 of the horizon volume. Also, as our main 
focus in this work is an order of magnitude estimate of iso-curvature fluctuations and since 
underline physical picture will remain same even for wiggly string, hence from now onwards, 
we only consider the case of straight string and the results will be quoted only for straight 
string case. 

We now study the density evolution of axions which may get trapped initially inside the 
wake. For that, let us hrst recall the effect of expansion of the universe on the dynamics of 
the phase transition. For that we mostly follow the work of Fuller [2] (see also, [6,7]) who 
have studied the density evolution of baryons inside the shrinking QGP bubbles. If R{t) is 
the scale factor of Robertson-Walker metric, then Einstein’s equations give [2], 

II - 

where p is the average energy density of the mixed phase and given by, 

P = M + (1 - f)ph (21) 

The energy density, pressure {pq,pq) in the QGP phase and in the hadronic phase(p/i, p^) 
are given by. 


Pq = gqaT^ + B, Pq = - B (22) 

Ph = PhaT^, Ph = (23) 

Here Pg ~ 51 and ~ 17 are the degrees of freedom relevant for the two phases 
respectively (taking two massless quark flavors in the QGP phase, and counting other light 
particles) [2] and a = ^. At the transition temperature we have Pq = ph which relates 
and the bag constant B as, B = ^aT^{gq — ph)- We dehne x = Pq/ph to be the ratio of 
degrees of freedom between the two phases. With this, Eq.(20) can be written as, 

+ ( 24 ) 

Now, conservation of the energy-momentum tensor gives. 


11 



( 25 ) 


During the transition, T and p are approximately constant. Therefore, Eq.(25) can be 
rewritten as. 


R{t) ^ f{x - 1) 
R(t) 3/(a; — 1) + 3 


(26) 


These two equations, Eq.(24) and Eq.(26) have to be solved simultaneously to determine 
the evolution of scale factor and the change of /, the volume fraction of QGP region as 
the transition proceeds. The initial time (to) relevant for us is when the overdensity inside 
the wake has been formed. Taking 15 long string inside the horizon at that epoch, we will 
consider the initial volume as, 

n «(^)4. (27) 

where dH{= 2to) is the size of the horizon at time to- Note that we take the wake like 
overdense regions to be well formed at time to. Thus our representative volume as the 
transition proceeds is, V (t) = ho(-^)^, Rq being the scale factor at to. Taking center of the 
wake as the origin and considering motion of the interface along direction, we can write 


= 2A{t)z{t). 


(28) 


Where, A{t) is the area of each sheet. Assuming the sheet extending across the horizon, we 
get the area as a function of time as, A(t) ~ E(t)^/^. Putting this value we get the evolution 
of the thickness z(t) as. 


z{t) 


m 

2 


° -Rn ' 


(29) 


Note that we are approximating the wake as bounded by two parallel sheets separated 
by average thickness of the wedge. Now let us determine the effect of interface motion on 
axion momentum distribution and subsequently on the evolution of number density of axions 
inside the wake. Total number of axions which will remain initially inside the wake can be 
calculated by integrating Eq.(lO). Thus, if N(to) be the total number of axions trapped 
initially inside the wake, then. 


r^max I^(t^ o 

iV(to) = Go / nkdki^)^ (30) 

^min -^0 

Where the last factor (^^)^ is due to the decrease of axion density from time i to to causes 
by the expansion of the universe (n^ as given in Eq.(lO) gives the number density of axions 
at time t). kmin and kmax are minimum and maximum comoving momentum respectively. 
Since, the minimum comoving momentum which can be trapped initially inside the wake 
cannot be smaller than (zq is the initial average thickness of the wake) we will take 

kmin as k or 2'kRo/zo whichever is larger. To determine that, we take the ratio. 
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( 31 ) 


2tiRq 


Rjf) ^0 _ / ^ \l/2 ^0 
i 2 tiRq ito 271 


Taking t from Eq.(9), to as ~ 10“^sec and the initial thickness zo ~ 1km or so, we get 
the ratio of 0(1). So, we take the minimum comoving momentum kmin as k. The maximum 
comoving momentum k^ax of the axions which will remain inside the wake can be set as 
discussed below. 


The momentum of the axions will keep increasing with the collapsing interfaces [9]. 
Thus, if k{t) be the comoving momentum of an axion whose initial momentum is ki, then 
the momentum at any time t, when thickness of the wake becomes z{t) can be obtained 
from the following relation. 


k(t) = ki 


zq R{t) 
z{t) Ro 


(32) 


The minimum physical momentum required for the axions to leave the interface is equal 
to Am (see, Eq.(l)). Therefore, the axions which have comoving momentum less than AmRo 
will remain inside the wake. Thus, we can put the upper limit of the integration in Eq.(30) 
as, k* or AmRo whichever is smaller. To determine that let us again calculate the ratio. 


1^* Rjt*) 

AmRo AmRo 



1/2 f 

Am 


(33) 


Putting the value of Am as obtained from Eq.(6) and Eq.(7) and the value of from 
Eq.(4) the ratio can be determined in terms of the formation scale of axionic strings as. 


k* 

AmRo 


- 7.9(i2-)-‘/^(--^; 

sec lO^^GeV 


~ 1.2 X 10^ 


Va 


W^GeV 


(34) 


Thus, the above ratio depends on the formation scale rja of the axionic string, which is 
constrained by the terrestrial and astrophysical experiments as well as from the cosmological 
considerations. The most stringent lower bound has been obtained from the SN 1987A [31] 
as Tja > 10^°GeE and the upper bound [15] (jja < iR^^GeV) comes from the consideration 
that the axions should not overdose the universe. If we take rja < lO^^GeE, then the above 
ratio is always less than unity. In this case we can take the upper limit of integration in 
Eq.(30) as k*. Having set the limits of integration we can now calculate total number of 
axions which will remain initially inside the wake and given by. 


iV(<„) = r„f|i-l](^)3 (35) 

Where, F = R{i)Grilln{riai)pc ^ :^r]lln{r]ai){ito)~^^‘^^. Now as the interfaces move 
towards each other momentum of each axion will be modihed according to Eq.(32). Let 
nRdk' be the modihed spectrum due to trapping of axions within the collapsing walls when 
thickness becomes z{t). Unless the momentum of each axion becomes AmR{t) the axions will 
remain inside the wake. So upto certain time ti, total number of axions N{t) (for to < t < ti) 
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will remain fixed inside the wake. Subsequently, number density will be increased due to 
decrease of volume fraction of QGP region / caused by interface motion. So, total number 
of axions at certain time t{< ti) can be written as. 


N{t) 




(to <t < ti) 


(36) 


Where, F' is new coefficient will be determined in terms of F. This coefficient F' es¬ 
sentially will take care of the change of the coefficient F due to shifting of momentum of 
each axion as is obtained from Eq.(32) and the change in QGP volume causes by interface 
motion. Determining this coefficient in terms F will be particularly useful in determining 
the evolution of axion density when axions will start leaking out the overdensity regions. 
k' = k^^^ = kc(t) (say) and k'^^^ have to be substituted by kc(t) and k*c(t), re¬ 
spectively. Substituting all these quantities and integrating Eq.(36) we get the total number 
of axions at any intermediate time between to to ti. Equating this number to the initial 
total number of axions at to as obtained from Eq.(35) one gets the coefficients F' as follows. 


VcFc^t) R(i) , 
fV{t) ' ’ 


( 37 ) 


Putting back F' in Eq.(36) and divided by the QGP volume fV(t) we get the number 
density of axions upto time ti and given by. 


pii) 


N(t) 

fV(t) 


y^F f(Q 3 1_ 

fV(tYR^’^k k* 


3 Vovl 
3271 fV{t)tl 


ln(T]at)[(ttof^‘^ 


(t*to)^^‘^[ 


(38) 


Now, as we mentioned earlier that the density written above is applicable upto time ti 
below which no axions will leave the QGP regions due to insufficient momentum required 
to escape the region. The thickness of the wake at that time ti can be calculated using the 
fact that when the momentum of axion exceeds the value AmR{t) they will leave the wake 
region, which is hrst happens when k*c(t) becomes equal to AmR{t). Using this equality 
we get the corresponding thickness at ti and determined by. 


Or, 


k'max = C(tl)k* 


zq R(F)R(t*) 
z(t) Rq t* 


AmR(ti) 


(39) 


z(ti) 

Zo 


Am 


(t*k 


)-l/2 ^ 7.8(iffi)-i/2( 


Va 


sec 


W^GeV 


Y ~ 1.2 X 10 


-3 


( 40 ) 


Since, the initial thickness Zo of the wake is ~ 1km, and to ~ 10~^sec, thickness of the 
wake at time ti comes out to be z{ti) ~ 2m for string formation scale r]a ~ lO^^GeU. As 
soon as the thickness decreases down to ~ 2m, the axions will start leaving the QGP region 
and finally all the axions will leave the wake and form a sheet like structure. To determine 
the evolution of density after ti, upper and lower limit of the integration of Eq.(36) should 
be replaced by AmR(t) and kc(t) respectively. Thus, the evolution of axion density within 
the wake after time ti is given by. 
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1 


(41) 


pif) 


VoFcjt) R{i) ^ 1_ 

fV{t) ^RitoY %{t) AmRit) 


3 Vovl 

3271 fV{t)tl 


ln{riai)[{ito) 


1/2 


Zq 

z(t)Am 


i,(«i <t<tj) 


( 42 ) 


Where, Eq.(37) and Eq.(38) have been used, t/ is the hnal time at which all the axions 
will leave the QGP region. The corresponding thickness can be obtained from the following 
expression, 

k'min = kc{tf) = = AmR{tf) (43) 

So, the ratio of the hnal thickness when no axion will remain inside to the initial thickness 
becomes. 


zjtf) 

Zo 


(44) 


Since, the time scale (fc) at the onset of QCD phase transition is of same order of magni¬ 
tude as our initial time R, which is of the order of ~ 10“^sec we get the ratio of hnal thickness 
to the initial thickness as ~ 7.7 x 10“^ for string formation scale rja ~ lO^^GeV. Therefore 
when the wake inside which axions were trapped reduced to a size of about ~ 0.1cm, the 
axions which were still remained inside the wake will all leave the wake and form a sheet like 
region. The release of the axions from the wake thus happens for the duration when wake 
thickness lies between 2 m and 0.1 cm. Below the thickness of 0.1 cm, no axion will remain 
inside the wakes. The axions thus left behind as the interfaces collapse will be concentrated 
in a sheet like regions. One can calculate the density prohle of these axions as follows. 


Let, Ppr{z) be the density of axions which is left behind at position z as the interfaces 
collapse. Then we can write down the following expression which relates total number of 
axions at position to the density which are left behind as follows. 


N{z) — N{z — dz) = Adzppr{z), 


( 45 ) 


where the time dependence of 2 : is given in Eq.(29). So We get, 

dN , , ^ ^ 

— = Appr{z). (46) 

Thus we hnally get the density prohle of axions as a function of thickness with collapsing 
interfaces as. 


Ppriz) 


^R{t) 



(47) 


However, here we should mention that in deriving the above equation we have considered 
that the axions will not disperse away immediately as it leaves the wake. Strictly speaking, 
the axions after leaving the wake will move with a velocity corresponding to the momentum 
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they have in the hadronic phase, as they leave the wake. Since axions leave as soon as 
their momentum in the QGP phase equals the mass difference Am (Eq.(l)), their kinetic 
energy may not be large in the hadronic phase. Hence they may not move very far from the 
wake (note, axions are very weakly interacting). In this case, the density will be somewhat 
decreased compare to the density as obtained from the Eq.(47). However the order of 
magnitude may not be changed much. So for simplicity we will consider the case where the 
axions will not be homogenized immediately after leaving the wake. In the following section 
we will discuss the results as is obtained from Eq.(36), Eq.(38) and Eq.(47). 


V. RESULT 

Eq.(24),(26) are solved numerically to get the evolution of scale factor and volume frac¬ 
tion / which is occupied by the QGP phase. The solution thus obtained is fed into Eq.(38), 
and Eq. (42) to get the evolution of axion density within the wake as the transition from 
QGP to hadron phase proceeds. Figl. shows the evolution of axion density inside the 
wake for two time zones. Time axis is given in the unit of psec, while density is in fm~^. 
As is shown in the figure, most of the time during which transition from QGP phase to 
hadronic phase {Attrans — 15/isec) happens, the axions will remain inside the wake. The 
reason is obviously due to dependence of axions density spectrum. The number of axions 
with minimum momentum will be more and they will leave the over density region last. In 
Fig.la, the density keeps increasing solely (for the time zone (to < t < ti)) due to decrease 
in volume of QGP phase. While for next time zone {ti < t <tf), there are two competitive 
processes. Decrease in volume causes to increase the axion density, while leaking out of 
axions will cause the density to decrease. Most of the time during which transition from 
QGP phase to hadron phase proceeds, hrst effect will be dominated over the latter, hence 
there will be net increase in density which is evident from the Figlb. Finally, almost all the 
axion will acquire the momentum needed to escape from the QGP regions is achieved and 
density sharply decreases down towards zero. The inset in Fig. lb shows the detailed prohle 
of decreasing part of the axion density plot. Fig2 shows the evolution of density (in unit of 
fm~^) with the thickness z((t) (in unit of meter). Fig2a, shows the plot upto thickness z{ti) 
until which density will keep on increasing due to the effect mentioned above. Fig2b, gives 
the evolution of density from thickness z{ti) to z{tf). For convenience, the inset in Fig2b. 
shows the expanded plot of the regions where density rapidly goes down towards zero. Fig.3 
shows the plot of total number of axions inside the wake as a function of thickness for the 
time zones, ti < t < tf. Here, volume of a wake is taken as the average thickness of the 
wake multiplying by the area A{t) (taking 15 long strings per horizon [24] and assuming 
the sheets extend across the horizon, area of each planar sheet is given as A{t) ~ as 
discussed earlier. The inset in Fig.3 shows the expanded plot within very small distance 
interval to illustrate the decrease of the total number of axions with the interface motions. 
Finally, the density prohle of axions which are left behind as the interface collapse which 
is obtained from Eq.(47) is shown in Fig4. Here, plot has been given for the relevant time 
interval (ti < t < tf). 

We can take average axion density of the order of ~ 10^ (in Fig.4), then we see that as 
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the interfaces bordering the wake is reduced to about ~ 0.1 cm, then almost all the axions 
will leave the wake (see Fig.3) and the density will be increased by a factor of ~ 10^. Since, 
the energy of the axions required to escape from the wake regions comes from the walls of 
the interfaces, these are kind of iso-curvature fluctuations. One can also calculate the total 
mass {Ma) of the axions which are concentrated within this planar sheet region (formed by 
left behind axions). Taking total number of axions Na of the order of ~ 10^^ (See Fig.3), 
the mass of the axionic sheet can be obtained from Eq.(7) and given by, 

M„ ~ Na ruah — 6.04 x 10“®eld(- )Na ~ 10^®(- )gm (48) 

rja Va 

Taking r]a — lO^^GeV, the upper bound on mass of the axionic sheet will be order of 
10^^ gm (~ at QCD epoch). Since, the fluctuations produced by axion trapping 

are of iso-curvature kind, they grow little [9,32] before radiation-matter equality era, teq. 
After teq, they grow in proportional to the scale factor. Hogan and Rees [32] have studied 
the evolution of iso-curvature fluctuations which were produced at the era of QCD phase 
transition. They have shown that iso-curvature perturbations produced by axions at QCD 
epoch can lead to formation of axionic ’minicluster’. The amplitude of the iso-curvature 
fluctuations produced in our model are very large and axions are concentrated within a very 
narrow sheet like regions of thickness of order 0.1 cm. This sheet will extend to a distance 
scale of order \/A{t). Where, A{t) ~ is the typical area of the wake as discussed 
above. This is about 2 km at QCD scale, which corresponds to comoving distance scale of 
order 10“^ Mpc today. One can study the evolution of such sheet like overdense region as 
it enters into the radiation-matter equality era. One can also study these overdensity at 
larger distance scales resulting from the large scale distribution of strings and their wakes. 
If these fluctuations survive until late stages, it will be interesting to study the effects of 
sheet like axionic clusters, especially whether they can have any effects on small scale CMBR 
anisotropies. 


VI. CONCLUSION 

In this paper, we have studied the iso-curvature fluctuations in axion density at QCD 
phase transition epoch due to the presence of density fluctuations produced by moving 
cosmic strings. We have considered the axions which are produced from the radiation of the 
axionic strings which are formed at some scale rja due to breaking U{1)pq symmetry. If the 
mass of the axions is relatively higher in the hadron phase compared to QGP phase then the 
axions may get trapped initially inside the wake-like overdensity regions. As the transition 
from hadronic phase to QGP phase proceeds with the motion of interfaces, these axions will 
acquire momentum due to collapsing interfaces and subsequently leave the wake. The axions 
thus left behind as the interfaces collapse may produce iso-curvature fluctuations. We have 
estimated the detailed prohle of the fluctuations. We have shown that, in our model the 
iso-curvature fluctuations in the axion density will be of order {^)axion ~ 10^ and they will 
be concentrated within a planar sheet like region of thickness few cm. This sheet can extend 
upto a distance scale of order 2 km at QCD scale. This length scale corresponds to comoving 
length scale of order 10“^ Mpc today. It will be interesting to study the implications of such 
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large fluctuations in the axion density especially on small scale CMBR fluctuations. Here, 
we should mention that, though we have discussed the density structure of wake for straight 
strings as well as for perturbed strings, results are quoted for straight string wakes only. For 
wiggly strings, dne to the lack of proper relativisitic analysis of shock formation, we have 
followed the analysis of ref. [22] and quote the resulting density flnctnations and opening 
angle of the wake (Eq.(16)) for wiggly string. The formation of shock in ref. [22] has been 
stndied for straight string moving through relativistic fluid. In their analysis, the flnid flow 
has been taken to be nniform. This assnmption may not hold for wiggly string. This is 
because, at the presence of wiggles, the fluid will experience the rapidly changing directions 
of the wiggles which causes acceleration of the flnid in different direction. So, in proper 
treatment of shock formation by wiggly string, one should take the non-uniformity natnre 
of the fluid in the analysis. We have also discussed the case where axion (being collisionless) 
directly (other than trapping mechanism) can give rise to density fluctuations. It turns ont 
that, the valne of these fluctnations are very small compared to the fluctuations produced 
from trapping of axions by cosmic string wakes. 
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FIG. 1. These figures show plot of evolution of axion density inside the wake as a function of 
time. The density p is given in unit of fm~^ and time t is in psec. Figures (a) and (b) show the 
density plot for two time zones to < t < ti and ti < t < tf respectively. 
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FIG. 2. These figures show plot of evolution of axion density p (in unit of fm~^) inside the 
wake as a function of thickness z (in meter) of the wake. Fig. (a) shows the plot upto thickness 
z{ti) and Fig.(b) is from z{ti) upto z{tf). 
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FIG. 3. This figure shows the decrease of total number of axions during the time interval 
ti < t < tf as a function of thickness 2 : (in meter). Small part of the apparently flat portion of the 
plot is expanded in the inset. 



22 




FIG. 

collapse. 


4. This figure shows the density profile of the axions which are left behind as the interfaces 
Density is given in unit of fm~^ and thickness z in meter. 
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